Bosonization in the two-channel Kondo model 
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The bosonization of the S = 1/2 anisotropic two-channel Kondo model is shown to yield two 
equivalent representations of the original problem. In a straight forward extension of the Emery- 
Kivelson approach, the interacting resonant level model previously derived by the Anderson- Yuval 
technique is obtained. In addition, however, a "(<t,t)" description is also found. The strong coupling 
fixed point of the (er,r) model was originally postulated to be related to the intermediate coupling 
fixed point of the two-channel Kondo model. The equivalence of the <t,t model to the two-channel 
Kondo model is formally established. A summary of what one may learn from a simple study of 
these different representations is also given. 



PACS numbers: 72.15.Qm, 72.15.Nj, 71.45.-d 

The physics of the isotropic two-channel Kondo prob- 
lem has been the focus of intense theoretical investigation 
since it was-iirst known to have an intermediate coupling 
fixed pointlll. The essential physics was understood at 
that time to be due to an overscreening of the impurity 
spin leaving a residual spin object to interact with the 
conduction electrons even in the strong coupling limit. 
This has since been clarified by the exact solution of the 
problem using the Bethe Ansatzotj. 

Considerable insight into this exact solution was re- 
cently provided by two apparently different approaches 
which both enable much of the physics of the problem 
to be extracted by elementary methods. In the first of 
these, Emery and Kivelsono found that a bosonized ver- 
sion of the model could be reduced to a resonant level 
model which becomes diagonal for a special value of per- 
pendicular coupling J z . The non- Fermi liquid thermo- 
dynamics and the Wilson ratio of 8/3 found near this 
pointll are consistent with the exact solution. In the sec- 
ond approach, Coleman et. al.il postulated a "cr + r" 
model, where an S — 1/2 impurity interacts with both 
the spin (cr) and charge (r) degrees of freedom of a sin- 
gle channel of conduction electrons. This model has a 
strong-coupling fixed point which, the authors argue, is 
equivalent to the intermediate coupling fixed point of the 
symmetric two-channel Kondo problem. A strong cou- 
pling expansion about this fixed point also gives a good 
account of the low energy physics of the exact solution 
and again reproduces the correct Wilson ratio. 

In analyzing the more general case where channel 
anisotropy in the couplings is allowed, a similar pattern 
of events has transpired. Andrei and JerezQ have ob- 
tained a solution for the thermodynamics of this model 
using the full machinery of the Bethe Ansatz. At the 
same time, an equivalent interacting resonant level model 
was obtained by Fabrizio et. alxi by the Anderson- Yuval 
technique and analyzed in the vicinity of the Toulouse 
point. The (s,t) model was also extended by Coleman 
and Schofieldcl to examine the fixed point of that model 
away from the symmetric point by a duality transforma- 
tion. Controversially, the authors argued that a Fermi 
liquid description of the strong coupling fixed point is 



not possible because of the Majorana structure of the 
Anderson model they derive at strong coupling. Up until 
now, the weak point in applying their argument to the 
two-channel Kondo model lies in the absence of the pre- 
cise link between this model and the (c,t) model which 
they treat. The need for such a connection has become 
all the more-pressing following recent papers by Zhang 
and Hewsor£3 which argue that the models have different 
fixed points on the basis of a comparison of the resistiv- 
ities of the two models. 

In this paper, I formally derive the (cr, r) description 
from the two-channel Kondo model using bosonization. 
The result is valid for both the channel symmetric and 
anisotropic cases. The resonant level model obtained by 
Fabrizio et. alS is also derived by this method (a result 
undoubtedly known by the authors of that paper). This 
gives a foundation for understanding how these simple 
pictures relate to the exact solution. In particular the 
explicit mapping between the (cr, r) model and the two- 
channel Kondo model shows why the "resistivity" calcu- 
lated by Zhang and HewsonE3 for the (cr, r) model has 
no correspondence with the two channel Kondo model. 
Recent work by Kotliar and Sitil has emphasized the im- 
portance of "Klein" factors in finding Toulouse points so, 
for completeness, I give a precise account of the bosoniza- 
tion technique including these factors. 

The two-channel Kondo problem describes how an im- 
purity spin (taken here to be S = 1/2) at the origin inter- 
acts with a conduction sea of spin-half electrons of two 
flavors (channels). Typically one restricts the interaction 
to being in the s-wave sector for each flavor reducing the 
problem to finding the radial wave functions of the con- 
duction sea. By mapping the incoming partial waves onto 
the positive x axis and out-going waves onto the negative 
x axis one can treat the problem as left-moving ID chiral 
electrons interacting with an impurity at the origintij 
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Summation is implicit over repeated channel indices v = 
1, 2 and spin configurations a =|, J,. S** and define the 
5 = 1/2 spinor at the origin. I have explicitly included a 
coupling to an external magnetic field in TIb and allowed 
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and 



asymmetry in both spin space [e.t 

between channels [e.g. jj 1 ^ ^ J±^]- 

I now proceed by bosonizing this model. This is a stan- 
dard procedure which I briefly review for chiral fermions 
moving a ring of length L with anti-periodic boundary 
conditions ip{L/2) — —ip(—L/2) (ensuring that kp = 
defines a unique ground state). Ultimately the limit 
L — > oo will be taken. As noted by Mattis and Liebt^l 
particle-hole excitations of a chiral ID model may be rep- 
resented by bosonic variables 
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These operators act on excitations above a vacuum state 
where all states below Uf are occupied and all those 
above are empty. Normal ordering '::' with respect to 
this vacuum is essential to avoid the infinities associ- 
ated with the infinite number of occupied states. Here 
q = 2T:n/L > 0. These bosons obey the correct commu- 
tation relations = 5 kq by virtue of the absence 
of a lower bound on the occupied momentum states. 
The electron creation operators are defined such that 

The idea behind bosonization is to use these bosons 
as an alternative to the fermion basis for describing ID 
models. The requirement that this new basis be complete 
means that, in addition to the bosons which describe 
particle- hole excitations, we also need ladder operators 
which create extra particles of a given spin and flavor in 
the Fermi seali-3, F i . These ladder operators (or "Klein 
factors" ) must also anticommute between different elec- 
tron flavors and spins to reflect the fermionic nature of 
the underlying particles created. They may also be de- 
fined to be unitary. There are many possible representa- 
tions for the ladder operators, and I chose to represent 



FJ by (- 



Here the Hermitian operator 9{ 



is conjugate to the number operator Nf. [6i,Nj] = £#y 
and [8i,8j] = 0. This formulation corresponds to a par- 
ticular ordering of the vacuum by the arbitrary assigning 
of an order to the spin/flavor degree of freedom i (e.g. 
1, | comes before 2, J,). By ensuring that the ladder op- 
erators commute with boson fields, Haldane was able to 
show that one may represent a Fermi field as 



if>l(x) = lim (— 1)^j< 



(7) 



where 

*i(x) = 0i + — + iJ2 

q>0 



(8) 



(a is a short wavelength cut-off which regularizes the 
momentum sums.) Together with Eq. 0, the following 
relations form a dictionary for bosonization 



[®(x)i, $j{y)] = -iSijirsgaix - y) 
d x $(x) = $'(ar) = 2tt : ^(x)ip(x) : , 
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Using this dictionary, I now express the Hamiltonian of 
Eq. [j] in terms of phase fields for each species of fermion. 
I choose {1 1,2 1,1 1,2 1} as my ordering convention 
for the Klein factors. Any ordering will work, but this 
one simplifies the computation as will be explained later. 
Following both Emery and Kivelson and Sengupta and 
Georges, the spin and charge degrees of freedom can be 
separated and both have symmetric and antisymmetric 
combinations across the two channels: 



c± = 7T [( $ U + $ li)± (*2 T +$2i) 



[(«1T - *u) ± (*2f - *2|)] 
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The number operators are defined analogously. 

In terms of this decoupling, the Hamiltonian may be 
written as 
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H.c.}, (15) 

S z & s _(0), (16) 
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From this Hamiltonian, I will first derive the interacting 
resonant level model obtained by Fabrizio et. aln. I will 
then show the connection to the (a, t) model. 

The resonant level model is obtained by the unitary 
transformation used by Emery and Kivelson 
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where a z is the 3rd Pauli spin matrix. The advantage 
of my original choice of ordering is that, since (— 1) s + 



2 



does not appear, this transformation does not affect the 
Klein factors. 

Using Eq. [ll] and the fact that, if [A, B] commutes with 
both A and B, then [e XA , B] = ALA, B}e XA , I find 

U*& a+ (x)U = $' s+ (x) - 2ttS z S(x) . 

With H' = WnU I have, H' = H Q and 
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Here I have written J z = (J^ + J^)/ 2 , <W Z = 
.1 - ')/2 and will similarly define J± and <5Jj_. 

To obtain the equivalent resonant level model I 
'refermionize' the problem by expressing the boson fields 
<£>a hi terms of fermion fields ip\. For example, 
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(21) 
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(Again the ordering of these fermions is arbitrary.) No- 
tice that since there are no terms in the Hamiltonian 
which change N c± , the operators (— 1) °± simply fix an 

arbitrary sign for Jj 1 " . Also note that for S = 1/2, I 

can define S + (— 1) s ~ = eft, where eft creates a fermion 
which anticommutes with ip\_ . (Maintaining the Klein 
factors is crucial to obtaining this relation.) 

Thus finally I may write the equivalent resonant level 
model for the two-channel anisotropic Kondo model as 



TL= ii'F 
J± 
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(24) 



The complex fermion d has been expressed in terms 
of its real and imaginary (i.e. Majorana) components 
d = (eft + d)/>/2 and b= (d} — d)/iV2. This emphasizes 
that the Majorana representation for the spin is perhaps 
the most natural since the conduction sea is seen to in- 
teract independently with the spin via its Majorana corria 
ponents. This is the model obtained by Fabrizio et. alu 
using the Anderson- Yuval technique — essentially identi- 
fying the perturbation expansions of the two modelsES. 
The charge degrees of freedom (A = c±) have completely 
decoupled and remain free — an observation that will be 
exploited in the second part of the paper. The Toulouse 



point corresponds to J z = 2nvp and 5J Z = where the 
model is non interacting and may be solved exactlyEj. 
Here the impurity spin completely decouples from the 
magnetic field so the local spin susceptibility is zero. The 
constraint 5J Z — means that in the Toulouse point 
only describes the region in the vicinity of the symmetric 
model. 

This model was analyzed in full in Ref. |s| for the re- 
gion around the Toulouse point. The two resonant level 
widths A b = Jj_/4nvpa and A Q = (5J±_) 2 /^nvp a, asso- 
ciated with the two independent Majorana modes, rep- 
resent the two relevant energy scales of the problem. 

A second equivalent description of the problem is ob- 
tained rather straightforwardly from the Hamiltonian 
expressed in terms of spin and charge degrees of free- 
dom Eq. 14-^. Refermionizing this Hamiltonian using 
Eq. B3, one obtains 



/°° d 1 f°° 

dx^l—^x + B S z + -J 



p s+ {x) : dx 



r (2) 



r(l) 



J; 



(2) 



+-|- : [ Ps+ (0) + Ps _ (0)] : : [p s+ (0) - p s _ (0)] : . 

The normal ordering of the density operators is an in- 
struction to subtract off the contribution coming from 
the vacuum state. The mean occupancy in the vacuum 
at x — is 1/2 for any fermion so I may replace : p(x) : 
by p(x) — 1/2. The (ct,t) description is apparent if I re- 
label subscripts s and c as channel indices 1 and 2, and 
subscripts + and — as f and |. This relabeling is shown 
pictorially in Fig. ^. I may then write the Hamiltonian, 
here for the rotationally invariant case, J± = J z , as 



1~t = Tto ^ 
+B 



j( 2 )<7 X (0) + J«Ti(0) 



■ S 



S+ cr\(x) + T\(x)dx 



where 



(25) 

(26) 
(27) 



This is precisely the model treated in Ref. || for the 
symmetric case [J^ 1 -* = J^) and in Ref. ^| in the gen- 
eral case [JW ^ j( 2 )]. The impurity is seen to interact 
now with the spin and charge degrees of freedom of a sin- 
gle chain leaving a second chain — representing what were 
previously charge excitations — decoupled from the spin 
(see Fig. [l]) . Since both the (cr, r) description and the in- 
teracting resonant level models can be derived from the 
original Kondo problem, they are both equivalent rep- 
resentations of the model. Thus, for example, the low 
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<t> + <[> 



■ft* 




(a) %+% (b) 
FIG. 1. A pictorial representation of the relationship be- 
tween the two-channel Kondo model (a) and the (a, r) model 

(b) . By expressing the boson fields $ in terms of spin and 
charge degrees of freedom (s±,c±) instead of the original fla- 
vor and spin decoupling (lj, j.,2f t |) I can refermionize the 
model to obtain the (a, r) formulation. 



energy spectrum of the (a, r) model plus the uncoupled 
chain will be identical to that of the two channel Kondo 
model. What then are the relative merits of these two 
equivalent models over the original formulation? 

By studying the non-interacting Toulouse point of the 
resonant level description a qualitative picture of the 
thermodynamics may be obtained. One sees a two-stage 
quenching of the entropy of the impurity spin as the 
temperature is lowered through the two energy scales. 
At temperatures between these two energy scales the 
thermodynamic response is controlled by the fixed point 
of the channel symmetric case. For example, the local 
spin susceptibility is logarithmically divergent as a func- 
tion of temperature. At low temperatures the specific 
heat and susceptibility revert to their Fermi liquid forms. 
These results mirror those obtained from the exact solu- 
tion. However, the disadvantage of this approach is that 
to solve the model one is restricted to the region near 



J. 



(i) 



( 2) 

J]' > Ji and so one may fail to identify prop- 
erties which are universal functions of the energy scales. 
For example, the Wilson ratio is vanishingly small near 
the Toulouse points in contrast to the rotationally invari- 
ant case where the Wilson ratio spans the range (0, 2]. 

It is in understanding the low temperature fixed point 
behavior that the (a, t) description excels. Unlike the 
original formulation, the strong coupling fixed point of 
the (a, t) model is stable: over-screening of the impurity 
spin is prevented by the Pauli principle. The intermedi- 
ate coupling fixed point of the two-channel Kondo model 
has beenmcwed to infinity by the transformation to the 
(a, t) picture. One may then exploit a strong coupling ex- 
pansion to examine the physics near the fixed point. Just 
as in the single channel caseEZrlia, symmetries restrict the 
form of the resulting fixed point Hamiltonian. As well 
as enabling a direct calculation of the Wilson ratio, one 
also obtains the form of the low-energy Hamiltonian. The 
Majorana character of this Hamiltonian is the basis for 
the conjecture that the excitations are not compatible 
with weakly interacting fermion quasi-particlesH. 

Finally I address the issue of how one uses the fixed 
point Hamiltonian of the (a, r) model to compute the re- 
sistivity in the two-channel Kondo model. The two mod- 



els are equivalent as we have shown but the mapping 
between them is linear only in the bosonized language 
of Eq. [l3| Put simply, the electrons in the two-channel 
Kondo model are not the same electrons appearing in 
the (a, t) model. To calculate the resistivity of the two- 
channel Kondo model one must show how the electron 
Green's function transforms under this mapping. It does 
not transform into the electron Green's function in the 
(a, t) model as implicitly assumed in Ref. [lC]. The appro- 
priate correlator is expressible only in terms of the boson 
fields: for example 



£i,t( x i 
1 



2na 



(28) 



where A = {c±, s±}. Following Affleck and LudwiglH, the 
low temperature resistivity in the isotropic two-channel 
Kondo model is found by considering this correlator with 
X2 < < x\ in the fixed point Hamiltonian of Ref. ^. The 
leading temperature dependence comes from considera- 
tion of the leading irrelevant operator which affects this 
correlator at first order and, since it has dimension 1/2, 
leads to a T 1 / 2 correctioncll 

In conclusion, two methods have been proposed re- 
cently to give a simple account of the physics of the two- 
channel Kondo model both at and away from the channel 
symmetric point. These methods sidestep, at the expense 
of completeness, the full formalism of the Bethe Ansatz 
by considering models related to the original two-channel 
Kondo problem. In this paper these related models are 
formally derived from the original two-channel Kondo 
model via bosonization. In describing the bosonization 
technique, care is taken to include explicitly the Klein 
factors which preserve the anticommutation relations be- 
tween the fermions on different channels. As expected, 
the interacting resonant level model found previously by 
Fabrizio et. al. is found as a simple generalization of 
the Emery-Kivelson Hamiltonian. In addition, however, 
the (cr, t) model with its strong coupling fixed point is 
also obtained directly from the two-channel Kondo model 
by bosonization. This establishes the equivalence of this 
model to the original two channel Kondo model. The 
relative merits of these two approaches are discussed. 
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